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The problem of the electromagnetic energy-momentum tensor is among the oldest and the most controver- 
sial in macroscopic electrodynamics. In the center of the issue is a dispute about the Minkowski and the 
Abraham tensors for moving media. An overview of the current situation is presented. After putting the 
discussion into a general Lagrange-Noether framework, the Minkowski tensor is recovered as a canonical 
' • energy-momentum. It is shown that the balance equations of energy, momentum, and angular momentum 

are always satisfied for an open electromagnetic system despite the lack of the symmetry of the canonical 
tensor. On the other hand, although the Abraham tensor is not defined from first principles, one can formu- 
late a general symmetrization prescription provided a timelike vector is available. We analyze in detail the 
variational model of a relativistic ideal fluid with isotropic electric and magnetic properties interacting with 
the electromagnetic field. The relation between the Minkowski energy-momentum tensor, the canonical 
energy-momentum of the medium and the Abraham tensor is clarified. It is demonstrated that the Abraham 
energy-momentum is relevant when the 4-velocity of matter is the only covariant variable that enters the 
■ constitutive tensor. 

On 

O \ 1 Introduction 

OO 

One century ago, in 1908, Hermann Minkowski U gave a complete relativistically covariant formulation 
of classical Maxwell-Lorentzian electrodynamics. In particular, he demonstrated how the macroscopic 
physical processes in moving bodies and media can be derived from the knowledge of the physics of the 
matter at rest. An important issue in phenomenological macroscopic electrodynamics of moving media 
is the definition of the energy and momentum of the electromagnetic field in matter. Rather surprisingly, 
this topic has demonstrated a remarkable longevity, and the question of the electromagnetic energy and 
momentum in matter was not settled until now, despite some theoretical and experimental achievements. 

It is not our aim to give a complete history of this issue. After Minkowski proposed (i I his version of 
the energy-momentum tensor, Abraham [2 , 3, 4) entered the dispute with a different expression. The main 
formal difference between the Minkowski and the Abraham energy-momentum tensors is that the latter is 
symmetric while the former is not. Since then, the discussion was going on which of these two tensors is 
correct (along with some other attempts, among which it is worthwhile to mention the results of Einstein 
and Laub 0, of de Groot and Suttorp J6| and of Peierls QUI)- 

Of the early studies, the work of Dallenbach and Ishiwara ||9j [10] QT| deserves to be mentioned. The 
later contributions of von Laue, Beck, Haus 021[nj[T4][l5] and many others ||T^ rTTl |T8l |T9l [20l [2T| [23l 
E2|23|26]|27]|28) and |[29l[33[3T][32l[33][34l[35l[3^[^ brought more arguments and more 

confusion into the discussion. Eventually, it became clear that the study of the electromagnetic field alone 
is insufficient, and the idea was put forward that all possible tensors are correct, depending on the way how 
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field and matter are separated into two subsystems, see for example [ 42 ,43 , 44 45] and l46l l47l l48l . More 
recent discussions were devoted to the analysis of this separation question 11491 l50l IBTl l52l l53l l54l 1331 . 
along with (56] E3 |5j |5j |gOj E0 12) and ESEIIS]- 

From a theoretical point of view, the most important developments should perhaps be attributed to 
Schmutzer and Schroder ll66l l67l l68l l69l who stressed the relation of the Minkowski and the canonical 
energy-momentum tensor, to Poincelot l70l I7T1 l72l |73| for underlining the universality of the Lorentz 
force, and to Penfield, Haus, and Mikura IT741 |75l l76l 1771 for coming up with the first working models for 
polarizable and magnetizable matter in the framework of the Lagrange variational approach. In our paper, 
these points will be put into the center of the discussion, following our earlier studies ll78l |79l . Some 
relevant results were also presented in the recent interesting papers of Dereli et al. If80l [8T1 . 

It is clearly impossible to mention even briefly all the names and to discuss all the contributions to 
the issue of the electromagnetic energy and momentum in a short paper. Moreover, since we will mainly 
pay attention to the theoretical questions here, the experimental results in this area will be not be properly 
discussed. In order to get more information, the readers should consult the thorough reviews [6 82 83 84 
[85j|86l|87l[83, see also the most recent one in |89l . 

Our basic notations and conventions follow the book 11901 . In particular, £o, Mo are the electric and the 
magnetic constant (earlier called vacuum permittivity and vacuum permeability). It is worthwhile to stress 
that we do not choose any special system of physical units. In other words, the presence of Eq and /io in our 
formulas does not mean that the international system SI is used, their values are different for SI, Lorentz- 
Heaviside, and rational systems, and they are fixed only after the specific unit system is chosen. More on 
physical dimensions and units can be found in |9T|. We do not discuss gravitational effects, the spacetime 
is flat. The Minkowski metric is g^ = diag(c 2 , —1,-1,-1). The totally antisymmetric Levi-Civita tensor 
rjijki is defined by 770123 = y/— 9, with g = detc/y . Latin indices from the middle of the alphabet label the 
spacetime components, i, j, k, ■ ■ ■ = 0, 1, 2, 3, whereas those from the beginning of the alphabet refer to 
3-space: a, b, c, ■ ■ ■ = 1, 2, 3. 



2 Preliminaries: macroscopic electrodynamics 

Quite generally, the Lagrangian of the electromagnetic field in matter can be written as 

L e = - ^ x m FijFu. (1) 

Here = c^A, — djA^ is the electromagnetic field strength constructed from the potential A ; , and 
the tensor density x 1 ^ 1 describes the electric and magnetic properties of matter. This quantity may also 
depend on Fij, thereby taking into account possible nonlinear electromagnetic effects, but in the first 
approximation ^ is a function only of the matter variables that describe the state of the medium. We 
will restrict ourselves to linear electrodynamics. 

As usual, we define the electromagnetic excitation tensor density by 

BT e 

For the theory ([T]i we then find the linear constitutive relation [92 .93, 94l l95l 

H ij = \ X ijkl F kt . (3) 

Accordingly, we then have L c = — j H kl Fki- 

In the conventional way, the components of the tensors Fij and IP* describe the electric and mag- 
netic fields (E, B) and the electric and magnetic excitations (D, H) (other names of which are "electric 
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displacement" and "magnetic field intensity"): 
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(4) 



It is worthwhile to mention that it was Minkowski HI who introduced the notion of a bivector (which he 
called "Raum-Zeit-Vektor II tcr Art", i.e., a second rank tensor in the modern terminology) and wrote for 
the first time the Maxwell electrodynamics in the explicitly covariant four-dimensional tensor form. 

In vacuum, the electromagnetic action reads I = —\J \og lk g-* 1 F^j Fki ^J—gd A x, where Ao = \J £o/mo- 
As a result, in vacuum we find for the constitutive tensor density 



X 



ijkl 



A 01 



-g(g ik 9 jl -9 il 9 jk )- 



(5) 



Consequently, the linear constitutive relation (O in vacuum in Cartesian coordinates reads (note that 

yf^g = c, Aqc = (Uq -1 , and A /c = e ) 



H ij = l%iJkl Fkl=F ij/^ 

or, equivalently, in 3 -dimensional form 

1 



D = e E, 



H 



B. 



/'o 



(6) 



(7) 



In macroscopic electrodynamics!], a polarizable and magnetizable medium is characterized by the non- 

p p 

trivial polarization P and magnetization M. They arise from the bound charge and current densities (p, j) 
via the standard relations 



p=-V-P, j=P + VxM. 
The vacuum constitutive law (0 in a medium is changed ||97l into 

1 



D = e E + P, 



H 



B M. 



/'o 



(8) 



(9) 



Analogously to the electric and magnetic fields, the polarization and magnetization are not 3-vectors, but 
constitute the components of the 4-dimensional polarization tensor of second rank, 
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Correspondingly, the constitutive relation (O can be recast into 
1 



(10) 



/'o 



p*3 +M l \ 



p. p p 

and, introducing the 4-vector of the bound current J 1 = (p, j ), equation ([8) can be rewritten as 



J 



djM ij . 



(11) 



(12) 



See Hirst 1961 for the definition and a discussion of magnetization and polarization in the microscopic theory. 
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Combining (O and (O, we can put the constitutive relation in an alternative equivalent form: 

M ij = - C jkl F kl , C jkl := x tjM - X ijkl ■ (13) 

In general, the susceptibility tensor density £, i: > kl describes also the magnetoelectric effects when the 
medium is magnetized in an electric field and/or polarized in a magnetic field. 

3 Abrahamization of an energy-momentum tensor 

Contrary to a general belief (and perhaps somewhat surprisingly) there does not exist any derivation of the 
Abraham energy-momentum tensor from first principles. In this section, we give a formal definition of 
what can be called the "abrahamization" prescription that can be applied to an arbitrary energy-momentum 
tensor. 

More precisely, let us assume that we can associate with a given physical system a second rank tensor 
ifc\ the components of which are interpreted as follows: to = w is a density of energy of the system, 
to a = s a (with a = 1,2, 3) is the energy flux density ("Poynting vector"), t a a = —p a is the momentum 
density, and t a b is the stress tensor. We do not assume that this tensor is conserved, and the divergence 
ditk = fk describes, in general, the balance equation of the system. For k = 0, this yields the familiar 
energy balance equation w + V • s = f . Integrating over a three-dimensional volume V, we find in a usual 
way that the change in time of the total energy J v w combined with the energy flux through the boundary, 
J gv s, is equal to the total power production of the system J v fa. Similarly, for = a = l,2,3we recover 
the balance equation — p a + db t a b = f a of the momentum of the system under the action of the force f a . 
The 4-vector force density vanishes, f k = 0, for a closed physical system that does not interact with other 
systems. 

Let us raise the first index, i u = g %k t\J> . The resulting covariant tensor is not symmetric in general: 
£*j ^ £?\ \ n particular, this means that the momentum density is not equal to the energy flux density, p ^ 
s/c 2 . There is nothing unphysical about this fact which may be explained, for example, by the nontrivial 
intrinsic angular momentum (spin) of the elements of the system [98 ). A well known symmetrization 
prescription of Belinfante and Rosenfeld 199"1|1001I101| allows one to construct from the original t k and 

B . B.. B .. 

from the spin density tensor another energy-momentum tensor t k which is symmetric, t 1J = t Jl , and 

b . 
tk l - 

A different symmetrization procedure is possible for the case when a timelike vector field u l is available, 
with u 2 = uiu 1 > 0. Physically, u l can be interpreted as a 4-velocity of the system with respect to the 
inertial reference frame. With the help of this field, we can introduce the transversal and longitudinal 

projectors P % - = 5j — u l Uj/u 2 and p*. = u l Uj/u 2 . Let us denote the antisymmetric part of the original 
energy-momentum tensor by a y := = — P l )/2. Then we define 

t k l :=t k * -P{ aj l -P)a k i. (14) 
By construction, this object is symmetric, 

t ij = - ^ u^u k a j)k , (15) 

and has the following crucial property: For the system at rest, u l = Sq, the time-time and space-space 

components coincide with the original symmetrized ones, to = to = w and t ab = t^ ab \ whereas the 
off-diagonal (time-space and space-time) components read 

t a0 = t a0 - 2a a0 = t 0a , t 0a = t 0a . (16) 
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A A 

In other words, the energy flux density ("Poynting vector") remains the same as before, s a = to a — to a = 

s a , but the momentum density p a is replaced with P a = t a0 = t 0a = s a /c 2 . 

The Belinfante-Rosenfeld procedure and the ab minimization ( TBI of the energy-momentum tensor both 
produce symmetric tensors from a given original asymmetric energy-momentum t^ % . However, these two 
symmetrization schemes are distinct in the following significant point. For the Belinfante-Rosenfeld ap- 

B . 

proach, the balance equation remains untouched because the divergence ditk 1 = ditk 1 is preserved. In 

A . A A . 

contrast, the divergence ditk 1 7^ d^k ''■ Let us find the difference of the two force densities, = ditk 1 
and f k = ditk' 1 . From ( [Pil l, 

h = fk - di{h<) - di{P)a^). (17) 

A 

In the rest frame, u % = <5q, we find fo = fo and 

f =f+ !( p -J0 +vxa - (i8) 

We use the boldface notation for the spatial 3-dimensional objects. In particular, f = {f a }, P = {Pa}, 
s = {s a }, whereas the antisymmetric part of the stress tensor gives rise to a = {\e a b c cL bc }- 

As such, the definition ( TBI looks rather artificial. To the best of my knowledge, there does not exist any 
derivation of ( TT~4-b from first principles, such as from the Lagrange-Noether machinery, for example. One 
merely demands the symmetry of the energy-momentum tensor under the condition that the energy density 

A . 

and the energy flux density in the rest frame remain the same. The form of the resulting tk x is then fixed 
by < TT~4T > which leads to the additional "Abraham force" terms (TT~8T > in the balance equations. 



4 Canonical energy-momentum for open systems 

It is sometimes claimed that the symmetry is a fundamental property of the energy-momentum tensor that 
is related to the conservation of the angular momentum of the system. In order to clarify this point, let 
us recall the relevant facts of the Lagrange and Noether theory. Suppose we have a system of the fields 
denoted collectively by $" 4 . Its dynamics is described by the action integral / = J Ld 4 x, where the 
Lagrangian (density) is a function of the fields and their derivatives, L = L{<^ A 1 di& ), and x l = (t, x a ) 
are the (local) spacetime coordinates. 

When the action / is invariant under the spacetime translations, one finds the conservation lavv0 

* E * i = -^pr** A - (19) 

Here the canonical energy-momentum tensor is defined by 
BT 

* i = m*r) 8 >* A - L *> m 

and we use standard notation for the the variational derivative, 

SL dL a ( dL \ 

- 9r T^r-TTT ■ (21) 



^ We omit the derivations of these results which are well known since the work of Noether [1021 ; for the detailed discussion see 
Sec. 19 of 11031 . for example. 
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Similarly, assuming that action is invariant under the Lorentz group, that acts on the coordinates and 
fields, we find 

Z w +d i S jk i = --^ I (p jk )£$ B . (22) 

Here (p-'i)^ are the Lorentz generators for the fields $ A , and the spin current density is introduced by 
FIT 

s ^-m^) {Plk)A ^ B - (23) 

Recall that the Lorentz generators are skew-symmetric, {pjk)^ = — (pfej)s- Hence the spin is skew- 
symmetric too, Sjk 1 = — Skj 1 '. 

When the physical system is closed in the sense that it does not interact with other systems, then its 
state is completely described by the variables <& A and their dynamics in time and space. The latter is 
determined by the field equations, 5L/5<fr A = 0. Using this in the conservation laws ( fT9] l and ( f22b . we 
obtain di £// = and Sum + <9; Sjk 1 = 0. Accordingly, we find that the canonical energy-momentum 
tensor is symmetric when the spin is trivial, Sjk 1 = 0. Otherwise, we can use the Belinfante-Rosenfeld |99l 
UOOIUOll procedure to derive the symmetric tensor Ok 1 = H^—dj + Sfc- 78 — Sk lJ ). By construction, 

we have di Uk = di £&% and err ifc j = 0. 

However, the above is true for closed systems only. For open systems that interact with other systems, 
we cannot put 5L/5$> A = and we have to keep the corresponding terms on the right hand sides of 
$1% and d22l . These terms describe forces and torques which result from the interaction of the systems. 
Even when the spin is absent, the canonical energy-momentum of an open physical system is, in general, 
necessarily asymmetric in order to maintain the balance of the nontrivial torque present on the right-hand 
side of ( f22l . We will see below how this works for the electromagnetic field interacting with polarizable 
and magnetizable matter. 



5 Energy-momentum tensors in electrodynamics 

After the preparations done in the previous three sections, we are now in a position to start the main 
discussion. Here we consider the construction of the canonical and symmetric energy-momentum tensors 
in the macroscopic electrodynamics of media. The corresponding Lagrangian approach was formulated in 
Sec. [2] Our starting point is the action I c = J Lfd^x. 

Without loosing generality, we can choose the collective field of the system as & A = (A*, The 
set of all fields is thus naturally divided into the two sectors, an electromagnetic and a material one. It is 
not necessary to specify how the tensor density x depends on the more fundamental material variables 
(we can always do this at a later stage), but instead it is convenient to treat x % itself as a generalized 
material variable. 



5 . 1 Canonical energy-momentum tensor 

The Lagrangian (fl} contains only derivatives of the electromagnetic potential A; but not of x » an d thus 
we easily derive from ( f2Qb and ([TJ the canonical energy-momentum tensor of the system: 

E fe 1 = - F k] - L c 5\ = - H ij F k3 + -W l F jt 5\ . (24) 

This is the well known Minkowski energy -momentum. Since the system is obviously open, the canonical 
tensor is not conserved. This becomes clear when we calculate the variational derivatives that enter the 
right-hand side of (1191 . Indeed, we find explicitly 

(25) 
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Both variational derivatives are clearly nontrivial. We should take into account that this (electromagnetic 
field) system interacts with the matter, and the total Lagrangian is thus the sum L tot = L e + L m . After 
introducing the electric current J 1 = SL m /5Ai, one then derives the Maxwell field equation 

r rtot 

— — = A' + J 1 = -BjH** +J l = 0. (26) 

This shows that A' vanishes only in the absence of electric charge and current densities of matter, but in 
general A 1 = — J 1 ^ 0. Substituting ( 1251 ) into ( fT9l ), we derive the energy-momentum balance equation for 
linear Maxwellian electrodynamics as 



diZk 1 = F kl J l + X k , X k = - F. l3 F mn d kX l ° mn - (27) 



The first term on the right-hand side is the familiar Lorentz force, whereas the second term requires some 
analysis. At first, we notice that 

A/. \ (FijdkH* - H ij d k Fij) = i (Fijd k M ij - M"ih F, i ) , (28) 
where we used ( fTTT ). Now, after some straightforward algebra, we can bring this expression into the form 

X k = -d i h i -F ki d j M ij , (29) 
after introducing the polarizational energy-momentum tensor 

L l := M l] F k] - - { M !l l- ,i 61 (30) 
The final step is to recall ( TTZl i and to substitute ( f29b into d27l i. After rearranging the terms, the result reads: 



5V + E fc *) = F ki ( J* + . (31) 



The form of this equation suggests to define the total electromagnetic energy-momentum tensor and the 
total electric current as the sums 

tot . e . p tot . p . 

Efe'-SV + lV, J l :=r + J\ (32) 

p 

The physical interpretation is clear: J % is the current density of the free charges, and J 1 is the polarizational 
current density of the bound charges. The corresponding energy and momentum, associated with the free 

c p 

and bound charges, are described by T, k l and T, k \ respectively. As a result, the balance equation ( Bit is 
recast into 

tot tot. 

^Efe* =F ki J\ (33) 

with the Lorentz force on the right-hand side that acts on all types of charges (free and bound) present in 
the medium. 

Combining (|24| | and (t3Qb . and using the constitutive relation (fTTT i, we derive the explicit form of the 
total energy-momentum: 

sV = (~H ij + M ij )F kj + ^(H jl - AP l )Fji 5i = ±-(- F ij F kj + ^F il F jt <SjJ . (34) 

This energy-momentum tensor was discussed by Poincelot ll70ll7Tll72ll73l and more recently by us |78| . 
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Quite remarkably, we discover that the form of the total electromagnetic energy-momentum tensor in 
media is precisely the same as in vacuum. It is worthwhile to stress that the final result (l33l ) is a direct 
consequence of the fact that the electromagnetic field is an open system ([TJ. The energy and momentum of 
the electromagnetic field are described by the canonical ( =Minkowski) tensor d24l) . but since the system is 
open, the balance equation d27l i contains nontrivial force terms that arise from the interaction of the field 
with matter. A careful evaluation of the force Xj~ then eventually brings the balance equation into the 
final form d33l l. This general result is the best what can be done without entering into the question about 
the structure and dynamics of the medium. In the next section we will discuss a specific model of matter 
and demonstrate how the additional knowledge of the physical nature of a medium can provide a different 
computation of the force term X^. 



5.2 Balance of the angular momentum 

The canonical energy-momentum (1241) is not symmetric. However, this does not mean that there is a 
problem with the angular momentum conservation. We simply have to recall once again that the system 
is open, and hence the right-hand side of the angular momentum balance equation d22b does not vanish, 
because the variational derivatives d25T l are nontrivial. 

The Lorentz generators for the material variable x rea d 

= + ^€s n M + SZS^X^, + SX^Sply (35) 

Using (f25t we then immediately find 



XTC 1 

|Vi ) m ™P? , , v m " p i =_£>., F , 4vi np9 = F<u, iff-," 



(36) 



We used (0 here. On the other hand, the antisymmetric part of the canonical energy-momentum d24l ) is 
straightforwardly found to be 

e 

^[jk] = -H[k n Fj]n = F[k\n\Hj] n - (37) 

We thus demonstrated that the balance equation d22l of the angular momentum, 

erg 

V\ , = ( „ \mnpq m'n'p'q' 

[jk] ~ Sx mnpq ' n ' p,q ' 

is satisfied identically for the Minkowski canonical energy-momentum. In physical terms, the skew part of 
the energy-momentum tensor is perfectly balanced by the torque which is present on the right-hand side 
because of the interaction of the field with matter. 

There is a somewhat subtle point which we silently avoided in our previous discussion. However, for 
completeness we have to mention it now. Namely, it is well known that there exist two possible choices 
of the electromagnetic variable: (i) We can take the 1-form A as fundamental variable, and construct the 
corresponding Lagrangian as a function of the field strength 2-form F = dA. This is our choice which 
we assumed in the book [90] and also here. The potential 1-form A is an invariant object independent of 
local coordinates and frames, and it does not transform under Lorentz rotations. The canonical energy- 
momentum 3-form then turns out to be explicitly gauge-invariant, given in tensor language by the expres- 
sion (|24| |. while the spin density formally vanishes, (ii) The second choice is to take the components A; 
of a covector. They transform (with a usual tensorial law) under the Lorentz group, and the corresponding 
generators read (p^k)™ = — $n d'S. For such a choice of the field variable, the energy-momentum tensor is 

c c 

given by the gauge noninvariant expression which is obtained from d24l ) by the shift Y^k 1 — * Sfc 1 + AE^ 1 , 
with an additional term AE^ 4 = —H^djAk. Since (p- J / c )™A m = — SnAu, we find from ( l23l the spin 
tensor Sjt- % = H l yA k ]. Its divergence thus contributes to the left-hand side of d22l which now reads 

AE W + n.S,,/ = - H^d^An + diiH^Aq) = (d z W b )A k] . (39) 
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But at the same time, a new term appears on the right-hand side of 

- A* = AyA fc] = {diH%)A k] . (40) 

Accordingly, we find identically 

xre 

AE [jk] + d t S jk * == - — (p jk )™A m . (41) 

Thus we come to the same conclusion that the angular momentum balance equation is perfectly satisfied 
despite the lack of the symmetry of the canonical energy-momentum tensor. 



5.3 Abraham energy-momentum tensor 

In Sec. [3] we described a general procedure for constructing a symmetric tensor from an arbitrary energy- 
momentum. Now we can apply this scheme to the electromagnetic field, taking the canonical energy- 
momentum ( T24l i as an input. Although we have demonstrated above that the argument of the "violation of 
the angular momentum conservation" is unsubstantiated for open systems, it seems reasonable to perform 
a detailed comparison of the energy-momentum tensors available on the market. 

In this approach, we assume the existence of the timelike vector field u l , with u 2 = c 2 . Then substituting 
( [24l and ( |37| > into the definition (fT4l . we construct the Abraham tensor of the electromagnetic energy- 
momentum: 



Ek* = --{H^F kj +F^H kj ) + -W l F jl 8l 
1 



2 [«*«, ( /• ,/. // "' - //,,/• ' ; ! + u k u l (/•'"//,, - // "/• .,;); . (42) 

A A 

In the rest frame, u l = Sq, we find the Abraham energy flux density s a = £o a and the field momentum 

A A o 

Pa = -Sa°: 

a s E x H 

P = - = —J—. (43) 

c c z 

The Abraham force (fT8l l then finally reduces to its well known rest-frame expression 

f = f+||DxB-^j+ivx(DxE + HxB). (44) 

The last word in deciding which energy-momentum tensor is correct (and in which sense) belongs 
certainly to experiment. However, the theoretical foundation of the Minkowski tensor obviously appears 
to be far more solid than that of the Abraham tensor. The former arises as a canonical energy-momentum 
tensor in the Lagrange-Noether framework, whereas the latter one is not derived from first principles. 
Nevertheless, as we will demonstrate in the next section, the Abraham tensor does resurface after we 
specify the structure and the dynamics of matter. 



6 Variational model of matter 

Let us now consider the dynamics of the medium. We will model the matter as an ideal fluid, the elements 
of which are structureless particles (i.e., no spin or other internal degrees of freedom are present). Such a 
continuous medium (see I11041I1051I106II . for example, for the relevant earlier work on the relativistic ideal 
fluids) is characterized in the Eulerian approach by the fluid 4-velocity «\ the internal energy density p, the 
particle density v, the entropy density s, and the identity {Lin) coordinate X. Furthermore, we assume that 
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the motion of a fluid is such that the number of particles is constant and that the entropy and the identity of 
the elements is conserved. In mathematical terms this means that the following constraints are imposed on 
the variables: 

diivu 1 ) = 0, (45) 
u%s = 0, (46) 
u l d % X = 0. (47) 

Due to the conservation of the entropy only reversible processes are allowed. In a variational approach, 
these constraints are taken into account by means of Lagrange multipliers. The classical action of the fluid 
reads I m = j L m d A x, with the Lagrangian density 

L m = - p{v, s) + A {u l u, - c 2 ) - vu l diki + A 2 u l d lS + A 3 u l d l X. (48) 

The Lagrange multipliers Ai , A 2 , A3 impose the constraints J4~5b-(l47b on the dynamics of the fluid, whereas 
Ao accounts for the normalization condition for the 4-velocity 

//,.</ V C-. (49) 

For the description of the thermodynamical properties of the fluid, the usual thermodynamical law ("Gibbs 
relation") is used, 

Tds = d{p/v)+pd(l/v) 1 (50) 

where T is the temperature and p the pressure. From this we have 

dp 
ds 

dp p + p 



vT, (51) 

as 



dv v 

One can in fact treat the above equations as the definition the temperature and the pressure. 



(52) 



6. 1 Equations of motion 

Recalling the general formalism of Sec. [4] we describe the material system (ideal fluid) by a collective 
field as § A = (u l , v, s, X, Ao, Ai, A2, A3). The last six variables characterize only matter, i.e., they enter 
the matter Lagrangian L m but not the Lagrangian of the electromagnetic field L c . The latter depends, in 
general, on the velocity of the medium u l and on the particle density v. This is a manifestation of the fact 
that both systems (electromagnetic and material) are open. 

The variational derivatives SL 1u /5Ak = 0, K = 0, 1, 2, 3, with respect to the Lagrange multipliers 
yield the constraints (|43T>-(f47T>. j49l , whereas variation of L nl with respect to s and X yields, respectively, 

^- = fl i (A a « < )+ — = 0, (53) 

OS c 

— = d i (A 3U i )=0. (54) 

In order to derive d53l ), we used d5Tb . 

It remains to find the variations with respect to v and u 1 . The direct calculation (use d52l ) yields 

^- = - P -±l-u%M, (55) 
5L m 

— - = 2k u i -vd i k 1 + k 2 diS + k 3 d i X. (56) 
du l 
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We cannot put these equations equal to zero because the material system is open. Contracting d56b with u l , 
after using d55l l and the constraints (l46b and ( f47b , we find the Lagrange multiplier: 

2Ao *(-,-,- J£. + Jfg). ,57) 



<- 



8v Su 1 



Substituting this back into d56l l, we derive the following useful relation 

-vd i A 1 +A 2 d i s + A 3 d i X = ^ lp + p + v — J +P{ -^— . (58) 

Here, as usual, P 1 - = 5j — u l Uj/c 2 is the transversal projector. Substituting this into d48l . we see that 
"on-shell" (i.e., if the equations of motion are fulfilled) the Lagrangian of the medium satisfies 

r) T m 

L m =p + v— . (59) 
ov 

6.2 Canonical energy-momentum of matter 
Since L m depends only on the derivatives of s, X, and Ai, we have 

n rm 

——j- d k <$> A = u l (-i/^Ai + A 2 d k s + A 3 d k X) . (60) 

Using (1581 and j59l , we then straightforwardly construct, from the definition ( f20b . the canonical energy- 
momentum tensor of matter: 

S fe 4 = u*P k + (-Si + -^u k u l ) p cS . (61) 



c 



Here we denoted 



p 8L m u k uJ SL m 
V k = ~^u k -\- — r - 5-^— r, (62) 

c z ou K cr ow 

P cfi = P + v^—- (63) 

The physical interpretation of these quantities is clear. The 4-vector d62l is the relativistic 4-momentum 
density effectively carried by the elements of matter. The first term on the right-hand side is the usual 
kinetic momentum determined by the mass (energy density) of the particles, whereas the two next terms 
arise from the interaction of the medium with the electromagnetic field. The same applies to the second 
term of the effective pressure (l63l l which "corrects" the usual hydrodynamical pressure by the term arising 
due to the fact that the material system is open. For the closed system we would have to put 5L m /Su k = 
and SL m /Sv = 0, and then the energy-momentum (l6Tb would reduce to the standard expression of the 
ideal fluid. 

Since the material system is open, the divergence of the energy-momentum tensor is nontrivial. The 
energy-momentum balance equation ( fT9l ) now reads 

m . SL ni 5L m 
d i ^ = -^d k u-—d k u\ (64) 

As usual, the right-hand side describes the forces that the electromagnetic field exerts on the matter. 

Let us now inspect the angular momentum balance equation. As a first step, we notice that only 
the 4-velocity vector field u l of the material variable <& A = (V, v, s, X, Aq, Ai, A2, A3) has nontriv- 
ial transformation properties under the action of the Lorentz group. The corresponding generators read 
{p^k)™ = fin Sfrt- However, since dL m /d(djU t ) = 0, the spin density of the medium vanishes, Sj k l = 0. 

Copyright line will be provided by the publisher 



12 



Yuri N. Obukhov: Electromagnetic energy-momentum 



The canonical energy-momentum of matter doTT ) is not symmetric, 

6L m 

However, the right-hand side of d22b is now 



We thus verify that the angular momentum balance equation (f22b is satisfied 



(66) 



m SL m 

*W s - (67) 
identically for the asymmetric canonical energy-momentum doTT i of the medium. 

7 Coupled system of the electromagnetic field and matter 

We finally can complete the picture by combining the two pieces which we studied separately above: the 
electromagnetic system of Sec. [2] and the material system of Sec. [6] In order to do this, we need one 
important additional input. Namely, we need to specify how exactly the two systems interact. In the most 
general form, the information about this interaction is encoded in the constitutive tensor density x m 
the electromagnetic Lagrangian (Q]). 

The model of the fluid, which we studied in the previous section, obviously describes the polariz- 
able/magnetizable medium with isotropic electric and magnetic properties. Such matter is characterized 
by the two scalar quantities: the permittivity e and the permeability //. In order to take into account the 
possible electro striction and magnetostriction effects, we allow for these quantities to be the functions of 
the particle density, 

e = s{v), [i = n(y). (68) 

We thus assume that the medium is intrinsically isotropic, and the possible anisotropic effects may arise 
only from the nontrivial motion of the medium, i.e., they are induced by the velocity of the fluid. 



7.1 Constitutive relation and the electromagnetic Lagrangian 

The constitutive relation for an arbitrarily moving isotropic medium is well known. It is given by the 
famous Minkowski equations, see Hi, Sec.E.4.2, eqs.(E.4.28)-(E.4.29) and (E.4.25)-(E.4.26). The con- 
stitutive tensor for this case is given by the formula that is very close to ©, namely, 



X 



ijkl 



X 



V wopt l aopti/opt aopti/opt 



where A = yiio7/i/io an d the so-called optical metric was first introduced by Gordon 111071 : 

1 - n 2 



i/opt — y 



(69) 



(70) 



Here n = -Jep, is the refraction coefficient of the medium. 

We straightforwardly find .g pt = c/n, and thus the Lagrangian (Q]i of the electromagnetic field 
finally reads 



Jj „kl rp P _ _ 



FijF 1 ' + 2- 



1 



F lk u k F il ui 



(71) 
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The components of the electromagnetic field F^ are given with respect to the laboratory system. In order 
to get some further insight into the formulas above, we notice that the electric and magnetic fields in the 
comoving system (in which the medium is momentarily at rest) can be described by the 4-vectors 

S,:=F lk A B':=^F jk u, (72) 
The field strength is uniquely reconstructed from these vectors as 

Fij = ^ (SiUj - SjUi + c Vljk iu k B l ) . (73) 
The electromagnetic Lagrangian ( TtTT > looks remarkably simple in these variables: 

L^-UeeoS 2 -— ), (74) 
2 V Wo/ 

with the obvious abbreviations £ 2 = £i£ l and B 2 = B t B l . Both vectors are by construction orthogonal to 
the 4-velocity, £ iU l = and B 1 ^ = 0. Hence, £ 2 < and B 2 < 0. 

In Sec. [5] we treated the whole constitutive tensor density ^ as the material variable. Now, after 
specifying the model of the medium, we have something better. The fluid under consideration is described 
by the variables u l and v, and the constitutive tensor now becomes a known function x = ^ (u m , v), 
given by the equation (f69b . The corresponding variational derivatives of the electromagnetic Lagrangian 
with respect to the material variables are easily computed: 

SL C n 2 - 1 . , 

-j— = jF kl F u u h (75) 

ou l MMoc 

SL C 1 ( ^ de c2 1 d^i 2 



\ £ °1T £ + 27T B ■ ( 76 ) 

bv 2 \ av ^oM v v J 

We used (|7T| i to derive (|75T l, however, it is much simpler to use a different (equivalent) form of the La- 
grangian (l74l to obtain ( |76| |. 



7.2 Canonical energy-momentum tensor 

After all these preliminaries, we are now in a position to find the explicit form of the energy-momentum 
tensors and to analyse the corresponding balance equations. 
The electromagnetic excitation (O for the Lagrangian ( |7TT i is 



H' 



_ ik jl p _ _ 
ffopt 9 opt r kl — 



n 2 - 1 



F^uuv? - F^uku 1 



(77) 



Substituting this into d24l i, we find the explicit form of the canonical (Minkowski) energy-momentum of 
the electromagnetic field: 



5V 



1 

n 2 - 1 



^^FflSl 



1 



F kn F nl u lU l - F kl u l F m u n + -F nl u l F^ Uj S l k 



This tensor is not symmetric, and its antisymmetric part (l37b now reads explicitly 

n 2 -l 



■u u F k]l F u ui. 



(78) 



(79) 
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7.3 Balance equations of the energy-momentum and angular momentum 

As we know, the Minkowski tensor ( 1781 is not conserved, and the energy-momentum balance equation < f27T > 
contains a nontrivial force on the right-hand side, X k = ^ FijF mn d k x % ^ mn ■ Since x^ M = x^ kl v), 
we have 

a ijmn a ijmn 

d k X ijmn = —~ d k v + -J— 8 k u m . (80) 

ov ou m 

Then the force density can be easily calculated, 



X * = ^-(lx ijkl F ij F kl )d k u+ 1 ^-^-F ij Uj F im d k u r ' 
ov \8 / MMoc 

= - S -f dk v- 5 -^-d k u\ (81) 
ov du l 

where we used ([D and ( 1751 ). 

Furthermore, for the explicit dependence of the constitutive tensor d69l l on u % and for the Lorentz gen- 
erators (f^k)™ = Sn S'S of the 4-velocity vector field, we prove the identity 

XTC ire 

(oi.) mnp9 Y mnpq =--—(oi) m u n (82) 

g x mnpq [ - fJ 3 k >m'n'p'q'X- - Su m 3 ' 

Making use of the variational derivative (J75J, we have 

XTC 2 _ i 

-S^MnU^-j^uyF^u, (83) 
Comparing with ( |79l ), we thus verify that the balance equation of the angular momentum 

EM=-|^fo*)n«" (84) 
is identically satisfied in this case, as before. 

7.4 Total energy-momentum tensor of the coupled system 

Now, the crucial step is to recall that the total system L tot = L c + L m of the coupled electromagnetic field 
and medium is closed. Thus, 

6L C + SL m _ Q SL C + 5L m _ Q 
Sv Sv Su l Su 1 

By combining (fsTT > with d64b . we then finally express the force that acts on the electromagnetic field in 
terms of the material variables, 

m . 

X k = -d^ k \ (86) 

Substituting this into (T27l i. and rearranging the terms, we find the true conservation law of the total 
energy-momentum of the closed system (field + medium): 

a, (fj fc * + £ fc <) = 0. (87) 

Note that we assumed that the medium is neutral, hence J 1 = 0. 
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Using the equations of motion ( I85l l. we can find the explicit form of the canonical energy-momentum 
of medium, too. In particular, combining d85l l with d75l ) and d76*l i. we have 



P 



-u k 



de 



F il Ul -^F jn u n F^ Ul ] 



1 2 



dv ^oA* 2 9u 

Remarkably, the effective pressure describes the electro- and magnetostriction effects. Substituting 
into d6Tl l. we obtain the final expression 



(88) 
(89) 



-6% + ^u k u^j p aS 



flflQC" 



(90) 



The total canonical energy-momentum tensor is the sum of the electromagnetic (Minkowski) tensor d78l 
and of the energy-momentum of the fluid d90l l. After some simple algebra we find 



u k u l + -5\ + -u k u l p on + 



n 2 - 1 



-F u u l F m u n 



F 3n u n F> l Ul + ±F nl u l F n i Uj 5i 



(91) 



It is worthwhile to note that whereas the Minkowski electromagnetic tensor and the energy-momentum 
of the fluid are both asymmetric, the total canonical energy-momentum is explicitly symmetric. This is in 
agreement with the fact that the total system is closed. Indeed, combining the balance equations (loTt and 
( T84T > of the angular momenta of the electromagnetic and material systems, we find 



\ Su m Su 

To derive the last equality, we used the field equations 



(92) 



7 . 5 Abraham energy-momentum tensor 
For completeness, let us compute the Abraham energy-momentum. In accordance with the definition < TT~4b . 

e c 

we find the projections of the antisymmetric part aj k = £uh of the Minkowski tensor d79l : 



r k u j 



P) a ^ 



■ u l F k 



U k U J 



F in F ui. 



(93) 



The resulting electromagnetic Abraham energy-momentum tensor reads 



A . 

Efe* 



5V - PiV - P\ a k i 



n 2 - 1 



-F kl u l F in u, 



1 

u k u l 



F ij Fkj + i F n Fjl si 

F jn u n F^ Ul + ]-F nl u l F^u 8l 



(94) 



The same expression is obtained, after some algebra, if we insert the electromagnetic excitation tensor ( |77| | 
directly in the general definition d42l . 

Comparing the formula (|94l with j9TT l, we observe that all the energy-momentum tensors are related via 
the interesting equation 



u k u' + I -61 + —u k u l I p cu + E fe 



„cff 



(95) 
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So to say, the Abraham tensor "absorbed" all the terms which explicitly contained the electromagnetic 
field, except for the electrostriction and magnetostriction terms, both from the Minkowski canonical tensor 
of the electromagnetic field and from the canonical energy-momentum of the matter. 

8 Why Abraham? 

The fact that the Abraham tensor resurfaced at the end of our derivations requires an additional analysis. 
Whereas the Minkowski tensor has a solid standing as a canonical energy-momentum deeply rooted in the 
Lagrange-Noether formalism, the Abraham tensor appears a rather artificial construct, as we saw in Sec. [3] 
So, why is it recovered in the equation (|95ll? Is this a coincidence, a specific feature of the models that 
we used for the description of the electromagnetic and material systems, or something more fundamental? 
Here we demonstrate that the latter is true, in a certain sense. 

8.1 A useful mathematical fact 

At first, we prove the following technical point. Let Qy = — ctji be an antisymmetric tensor and /3k a 
4-vector. Then cty satisfies the algebraic equation 

if and only if this tensor is constructed from /3k and the velocity u % as 

otij = /3[i uj]. (97) 
The proof is straightforward. Let us assume that d97l ) is true. Then 

H< = {si-^f)l(0 j u i -P i u j )=^^/3 j , (98) 
iW = ^5T \ i^u j - (3 J u k ) = y p{ f3j. (99) 



Hence, ([96) is fulfilled. 

Conversely, suppose d96l l is true. Contracting this equation with u i , we find 

I'i UiOj* + UiP) a k i = 2 u t a k l = c 2 I'i (3 j =^ v?uf = - °- p\ $ . (100) 
On the other hand, from ( |96l l we have, making use of ( 1 100b : 

I'i otj* = -P) ^ ] +u?PiPj = \ ^Plh- (101) 
This finally yields, again using ( 1 100b . 

= (H + ^ - \ {f Pi ft - ft ) - i (ft* - -ft) ■ (102) 
Thus, the statement is proved. 

8.2 Crucial role of velocity 

We now extend the electrodynamical model of Sec. l7.ll bv generalizing the constitutive law j69l . Namely, 
we assume that the constitutive tensor is not a particular function d69l of the 4-velocity u l and of the 
particle density v, but it rather is some general function 

X ijkl =X ijkl (u m ,P,ilj n ). (103) 
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Here the collective variable ip n denotes all possible parameters that describe the state of the polariz- 
able/magnetizable medium. In the macroscopic approach, ip includes the most general matrices of di- 
electric, magnetic, and magnetoelectric susceptibilities. An important assumption is that these parameters 
are not changed under the Lorentz transformations relating different reference systems. 

More precisely, we thus assume that the 4-velocity of the medium u 1 is the only (Lorentz) covariant 
object that enters the electrodynamical constitutive relation ( 1103b , and the other material variables (v, ip ) 
are all scalars under the Lorentz transformations. Technically, this is realized by assuming that the parame- 
ters ifr 1 take their genuine (or intrinsic) values in the comoving reference system with respect to which the 
medium is momentarily at rest. In other reference systems, these intrinsic values of ip n remain the same, 
whereas in every reference system the form of the constitutive matrices is determined only by the motion 
of the medium, i.e., in technical terms, by the velocity u % . 

This assumption has the far-reaching consequences. Since the velocity u % is the only covariant argument 
(with the Lorentz generators (p-'fc)™ = 5n of the constitutive tensor x* 3 ' fc ' v, ip Q ), we have 

A T c , A T e r)-v mn Pi A T e A T c 

U±J I „ \mnpq ^m'n'p'q' _ L '■ LJ ~A ( „ \m „,n _ ~~_ ( „ \ro„,n _ U±J „, nf\A\ 

Sx mn Pq (Pjfc)m'nW* ~ du m UWn " ~J^^ k ^ U ~ fo[k U o\ ' ( 104) 

[Hint: denoting u jk = -lu^ the parameters of the Lorentz group, i^™ M = u jk (p jk )2Z? P > q >X m ' n ' p ' q ' 
gives the infinitesimal Lorentz transformation; using for differentiation the chain rule completes the proof.] 
As a result, we notice that the angular momentum balance equation ( T3"8l l is of the form d97l ) where 

Oij=V m , Pk = j^- (105) 
Accordingly, combining d96l l with the definition of the Abraham tensor ( fT~4T >, we find the relation 

e . A . . X ■ 5L e 

Vk l = Y>k l + u l p{—. (106) 

ou 3 

Adding the canonical energy-momentum doTb . d62b . we then finally come to the equation d95l ). 

This analysis demonstrates that the relation d95l l between the canonical energy-momentum tensors and 
the Abraham energy-momentum is not occasional. It is valid not only for the isotropic constitutive relation 
d69l ) of the specific model which we studied in Sec. 17.11 but for the general constitutive relation d 1 03t as 
well. Sec. l8.1l provides an explanation why the abrahamization prescription dT4b is successful, after all. 

The key to the Abraham tensor is in the 4-velocity of the medium: when u % is the only (Lorentz) 
covariant variable that enters the constitutive relation (and only then), the projections of the Minkowski 
tensor satisfy d96b that yields ( 1106b . This explains why the Abraham tensor turns out to be relevant for 
the discussion of the energy and momentum of moving media, especially of the isotropic ones. However, 
the angular momentum balance equation no longer has the form d97l ) if the constitutive tensor depends on 
additional covariant material variables, and the final nice result ( 11061 ) is invalid then. The relevance of the 
Abraham construction is thus limited. 



9 Discussion and conclusion 

The complete century-long history of the discussion of the electromagnetic energy and momentum in 
moving media is still to be written. In this short paper, we presented a consistent viewpoint on this problem 
within the framework of the Lagrange-Noether approach. The canonical energy-momentum tensor appears 
then as a fundamental structure, together with the relevant balance equations of energy, momentum, and 
angular momentum. 

By a detailed derivation, we explicitly demonstrate that the angular momentum balance equation is 
always perfectly satisfied for the canonical tensor of the electromagnetic system when one takes care of 
the fact that it is an open system. Even in absence of the intrinsic (spin) degrees of freedom, the canonical 
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energy-momentum tensor of an open system does not need to be symmetric, contrary to a widely spread 
misunderstanding. 

The Minkowski energy-momentum tensor d24b arises as a canonical structure in macroscopic Max- 
wellian electrodynamics, which stresses its fundamental physical nature. The balance equation of an open 
electromagnetic system (l27l i contains nontrivial force terms. Without specifying a model for the medium, 
we demonstrate that the Xk force d28l > is related to the bound charges and currents, thus giving rise to the 
balance equation d33l of the total canonical energy-momentum tensor d34l >. Its validity for the analysis of 
the crucial experiments was demonstrated in ll78l . 

Alternatively, the Xk force (l28l l can be reconstructed in terms of the material variables when a model 
of the medium is specified. We show how this works in the framework of the variational approach to the 
relativistic ideal fluid. Here again the canonical energy-momentum doTb appears as a fundamental object. 
The interaction of matter with the electromagnetic field induces an additional field-dependent term in the 
material momentum d62l i superimposed with the usual hydrodynamical structures. The pressure d63l picks 
up the electrostriction and magnetostriction contributions. 

In contrast to the canonical structures, the Abraham energy-momentum tensor does not arise from first 
principles. We give a formal prescription that allows one to construct a symmetric tensor ( TT~4T > from any 
given energy-momentum, provided a timelike vector field u l is available. Unlike the other symmetrization 
procedures (such as the Belinfante-Rosenfeld one, for example), the abrahamization does not preserve the 
balance equations, it rather introduces the additional forces in (TTTt . dT~8b . When applied to the Minkowski 
tensor, this prescription yields the famous Abraham electromagnetic energy-momentum d42l) . 

Nevertheless, despite its ad hoc definition, the Abraham tensor turns out to be relevant since it resurfaces 
at the end from the sum of the two canonical energy-momenta of the electromagnetic and material systems 
in the eq. ( |95l l. It thus appears as a certain "hybrid" construct that absorbs the explicit field-dependent 
terms from the canonical energy-momenta of the electromagnetic and material systems. This curious fact 
explains why the Abraham tensor was for such a long time seriously competing with the Minkowski tensor 
in the analyses of the numerous experiments ll88l . 

This result is not confined to the specific model with the constitutive relation d69l l. but holds in general 
for the class of constitutive tensor densities ( 11031 ) that depend arbitrarily on the material variables, provided 
the 4-velocity u l of the medium is the only Lorentz-covariant object. This apparently confirms the recent 
observations of Dereli, Gratus and Tucker ll80l [BP that the Abraham energy-momentum arises as the 
metrical energy-momentum tensor for media with a general constitutive relation. If we assume that the 
metric of spacetime is not the flat Minkowski one but is a function of coordinates, the metric energy- 
momentum tensor density is defined by the variational derivative 

a v :=2c ^ L ] +Lm \ (107) 

e m 

Indeed, a (long but straightforward) computation then shows that oV 1 = Ylk 1 + Y>k for the model (Qji, 
(f69t and (l48l l of the isotropic moving medium above. There is a difference, strictly speaking, between our 
results since neither the Lagrangian of the fluid nor the hydrodynamical part of the energy and momentum 
is discussed in [|80ll8Tl . 

It is worthwhile to stress that the model of the medium above is by no means a general one. Among 
its most serious limitations is that the elements of matter are assumed to be spinless. The extension of the 
model to the case of the nontrivial spin (along the lines similar to those of l46l l47l l48l ) can bring new 
insight into the subject, especially with regard to the effects of magnetism. 
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